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Clusters of spherical particles with isotropic attraction favour compact structures that maximise the number of energetically 
optimal nearest-neighbour interactions. In contrast, dipolar interactions lead to the formation of chains with a low coordination 
number. When both isotropic and dipolar interactions are present, the competition between them can lead to intricate knot, 
link and coil structures. Here, we investigate how these structures may self-organise and interconvert in clusters bound by 
the Stockmayer potential (Lennard- Jones plus point dipole). We map out the low-lying region of the energy landscape using 
disconnectivity graphs to follow how it evolves as the strength of the dipolar interactions increases. From comprehensive surveys 
of isomerisation pathways, we identify a number of rearrangement mechanisms that allow the topology of chain-like structures 
to interconvert. 
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1 Introduction 

The field of chemical topology : can be traced back to the first 
characterisation of interlocked organic rings and the introduc- 
tion of the term catenane by Wasserman in I960. 2 The two 
rings in a catenane are not chemically bonded, and the topol- 
ogy must therefore be specified to distinguish the structure 
from the unlinked rings, giving rise to the concept of topo- 
logical isomers. Wasserman's short report also points out that 
a single cyclic hydrocarbon of sufficient length may, in princi- 
ple, exist not only in the form of a simple loop but instead as 
a topologically distinct knot. 

Since this early work, there has been considerable progress 
both in the theoretical understanding 34 of "mechanically 
linked" cyclic molecules and in efficient methods for syn- 
thesising them. 5 Knotted topologies have even been discov- 
ered in organic synthesis even when originally unexpected. 6 
The prevalence and importance of knots in biological macro- 
molecules has also been recognised. Catenated and knotted 
forms of circular DNA naturally occur, for which topological 
isomers are enzymatically interconverted during replication, 
transcription, and recombination. 7 8 Knot-like motifs arise in 
the native structure of a surprisingly large number of proteins. 
Although proteins are unbranched chains — and therefore have 
a topology that is formally trivial — considerable insight into 
the role of knots can be gained by chemically joining the ends 
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of the protein to make a cyclic topology that contains a gen- 
uine topological knot. 9 

Cyclic structures can also emerge from the self-assembly of 
particles that interact via non-covalent interactions. In particu- 
lar, colloidal particles that carry a dipole moment tend to form 
chains that can close into loops, 10 and the possibility of knot- 
ted topologies in such systems has been pointed out. : l Highly 
intricate knotted topologies also arise in other types of col- 
loids in a way that is conceptually different — not from chains 
of particles themselves, but from lines of defects in cholesteric 
liquid crystals, where twisted nematic order is disrupte d by 
anchoring to spherical particles suspended in the solution JEEE1 
Once again, however, it is the anisotropic interactions between 
particles that leads to the existence of these lines and the rele- 
vance of topological considerations. 

Global optimisation calculations have shown that knots, 
links and coils are energetically favourable structures for 
clusters of Stockmayer particles, which possess a permanent 
dipole plus an isotropic soft core and attractive tail. 14 These 
intricate topologies arise from the competition between the 
dipole-dipole interactions, which promote chain formation, 
and the isotropic attraction, which favours compact, highly- 
coordinated structures. The energetically optimal compromise 
is often a continuous closed-loop chain, which is entwined 
with itself or with another loop in order to increase the number 
of nearest-neighbour contacts. 

Most of the literature concerning the Stockmayer model has 
addressed the structure and thermodynamics of the bulk fluid 
phase. In this context, too, the interplay between chains and 
compact structures is crucial. Van Leeuwen and Smit were 
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amongst the first to point out that gas-liquid phase separation 
is suppressed if the dipolar attraction dominates. 15 Ten Wolde 
et al. showed that droplets of Stockmayer particles appear 
from the vapour phase first by formation of chains and then 
by collapse of the chains into globules that nevertheless re- 
tain a chain-like structure. 16 The more general field of dipolar 
fluids, including the Stockmayer model, has been summarised 
in an informative review by Teixeira et al. 17 There has been 
relatively little work on finite clusters bound by the Stock- 
mayer potential, 18 though the model has now been recognised 
as having applications in self-assembly. 19 

The potential for a cluster of N Stockmayer particles takes 
the form 
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where T\j is the position vector of particle i with respect to 
particle j, £; is a unit vector along the dipole moment of par- 
ticle /, £ and a are the Lennard- Jones (LJ) units of energy 
and distance, respectively, and jl is a dimensionless parame- 
ter that determines the relative strength of the dipolar and LJ 
components of the potential. As fl increases and the dipolar 
contribution to the energy dominates, the Lennard- Jones well 
depth £ is no longer a convenient unit of energy. Here, we will 
report energies in units of the well-depth £* of the full Stock- 
mayer pair potential for parallel head-to-tail dipole vectors at 
the relevant value of fl. 

The structure of the clusters is tuned by the dipole moment 
strength ji. In the limit of a weak dipole (small fl), com- 
pact, icosahedral structures are favoured. 20 For strong dipoles 
the particles behave as dipolar soft spheres, producing chains, 
rings, and branched clusters, such as those formed by assem- 
blies of ferromagnetic nanoparticles.^It is in the intermediate 
regime, where competition between the isotropic and direc- 
tional contributions leads to frustration, and interesting topolo- 
gies with as many as ten irreducible crossings may arise. 14 

The current work builds on these findings concerning en- 
ergetically optimal structures by making a broader survey of 
the energy landscape for a selection of cases. The analysis 
presented here enables us to characterise the overall evolu- 
tion of the energy landscape between the Lennard- Jones and 
dipole-dominated limits, and to identify pathways and rear- 
rangement mechanisms by which different morphologies in- 
terconvert. Some of the key properties can be observed in 
the small St 13 cluster (where St^ denotes an assembly of N 
Stockmayer particles). The analysis is then extended to more 
complex rearrangements of knotted morphologies in the St2i 
and St38 clusters. 



2 Methods 

2.1 Optimisation 

A Stockmayer particle consists of a single interaction site with 
an isotropic Lennard- Jones component and an anisotropic 
point dipole. Since the dipole has cylindrical symmetry, each 
particle has five degrees of freedom: three translational and 
two rotational. In previous work 14 on the global optimisation 
of Stockmayer particles, the translational coordinates were 
represented in Cartesians and the orientation of the dipole was 
described using the polar angle and azimuthal angle (f) . 

The present work requires optimisation of first-order saddle 
points and the calculation of pathways, as well as characterisa- 
tion of minima on the potential energy surface (PES). For such 
applications it is more convenient to represent the orientation 
of the particles using the general angle-axis framework J 22 ! 23 1 
where the orientation of a rigid body is obtained by rotating a 
particle from a reference orientation about an axis specified by 
the rotation vector p through an angle given by its magnitude 
|p|. Although this representation introduces a redundant sixth 
degree of freedom for each particle in this case, it has the ad- 
vantage that the components of p can be treated as unbounded 
variables, unlike and (j). The redundant degrees of freedom 
contribute additional zero eigenvalues to the Hessian matrix; 
the corresponding eigenvectors can be obtained analytically^ 
and are used for projection while characterising the transition 
state as described later. 

In the case of Stockmayer particles, we arbitrarily took the 
reference orientation of the dipole along the laboratory-fixed 
z axis. The unit dipole vector pL i of a particle can then be 
obtained from the angle-axis variables p;. All angular first 
and second derivatives of V can then be written with respect 
to the components of p/. 

2.2 Calculating pathways 

A potential energy landscape can be characterised by mapping 
out its local minima, which correspond to locally stable inher- 
ent structures,^ and their connections via first-order saddle 
points (transition states). 

The survey of an energy landscape typically starts from a 
small selection of energetically low-lying minima obtained 
from a basin-hopping p 5 ! 26 ! global optimisation calculation. To 
find the pathways that connect pairs of these structures, 27 we 
employ the doubly-nudged elastic band method, 28 as imple- 
mented in the OP TIM packaged In this approach, a discre- 
tised interpolation is first established between the two min- 
ima in the full configuration space. The path is then allowed 
to relax on the PES while the states are kept approximately 
evenly spaced along the path by springs. Local maxima on the 
path are then tightly converged to transition states using hy- 
brid eigenvector- following, E3E2 and the connectivity is deter- 
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mined by approximate steepest-descent minimisations, adding 
any new minima thus found to the database of structures. It is 
common for new minima to be found during this procedure 
and so the first iteration of attempted connections rarely pro- 
duces a fully connected path between the minima that were 
originally specified. 

Missing sections on a path are approached in the same way. 
Candidate pairs of minima to be bridged are selected by re- 
peated application of the missing connection algorithm, 27 also 
implemented in the OP TIM package, 29 until a connected se- 
quence of minima and transition states (a discrete path) is 
found. To characterise complete paths, the minima are consid- 
ered to be nodes on a graph and the weight of the edge between 
two minima is set to zero if they are already directly connected 
by a single transition state in the database. If no connection is 
known, the weight is set to a function of the shortest Euclidean 
distance between the minima. 27 Hence, the algorithm identi- 
fies breaks in the path that are short in configuration space, 
making it more plausible that a connection will be found. The 
Euclidean distance between different Stockmayer clusters was 
determined from the positions of the particles only, ignoring 
the orientations of the dipole moments. 

2.3 Discrete Path Sampling 

Once an initial path has been established, the stationary point 
database must be expanded in order to obtain a more complete 
picture of the landscape. For small systems, a comprehensive 
survey of the landscape is possible by repeatedly attempting 
to connect all pairs of minima already in the database by the 
method described in Section [Z2| until no new minima or path- 
ways emerge. This approach was employed for the St 13 clus- 
ters. 

For clusters even just a few particles larger, it is not feasible 
to obtain exhaustive lists of stationary points. It is then most 
informative to obtain a good representation of the regions sur- 
rounding the lowest-lying minima of each structural family 
and a set of kinetically relevant pathways between them. Al- 
gorithms for connecting regions of configuration space in this 
way using discrete path sampling 33 35 (DPS) have been de- 
scribed in detail in previous work.ESEH Consider two regions 
on the landscape, A and B, which contain the minima we wish 
to connect. Minima in neither region that nevertheless appear 
on pathways between A and B minima belong to the interven- 
ing region, denoted I. Assuming Markovian dynamics, that 
local equilibrium is reached within the A and B sets, and that 
the steady- state approximation applies to minima in the I set, 
the steady- state rate constants &|| and kj^ can each be writ- 
ten as a sum over discrete paths. 33 In turn, the contribution of 
a given discrete path can be expressed in terms of the individ- 
ual rate constants k(j for each successive pair (/, j) of minima 
along the path. 



Since we are concerned with characterising the PES rather 
than obtaining quantitative values for the rate constants, we 
estimate the single- step rate constants ktj using a simplified 
version of harmonic transition state theory, where the ratio of 
normal mode frequencies at the transition state and the mini- 
mum is set to unity. This simplification has the advantage of 
making it unnecessary to consider the relationship between the 
vibrational and librational contributions to the normal modes 
of the Stockmayer particles. To proceed otherwise would re- 
quire specialisation of the model by choosing a particular mo- 
ment of inertia for the particles. With this simplification, the 
expression for the rate constant reduces to 



*y = -Yexp[-(v;t-v;-)/(tor)], 



where 0[ denotes the order of the point group for minimum /, 
and the superscript f is used to refer to the transition state. For 
the calculation of rate constants, we work at a temperature of 
k^Tfe* = 1/30, which is the regime probed by some recent 
experiments in which rings of magnetic dipolar colloids have 
been observed. 21 

Using the simplified rate constants for individual steps, the 
discrete path that makes the largest contribution to k^ A (the 
fastest path) can be identified. Attempts are then made to 
shortcut this path by directly connecting minima that lie on it 
but are separated by intervening minima. The choice of min- 
ima can be made on the basis of their separation along the path 
or in Euclidean space, or from the height of the barrier that 
separates them. These schemes guide the exploration of the 
PES towards kinetically relevant minima and transition states, 
and are here applied to the more complex landscapes of the 
St2i and St38 clusters. 

It is possible for artificial kinetic traps to arise as the 
database is generated — that is, for minima to be connected 
to the A and B regions by large barriers, where smaller bar- 
riers exist but have not yet been found, resulting in spuriously 
small rate constants. As the database grows, we periodically 
attempt to remove such traps by identifying pairs of minima 
based on the ratio of the potential energy barrier (from min- 
imum to transition state) to the potential energy difference 
between the two minima connected by the transition stated 
These pairs are then subjected to reconnection searches as de- 
scribed above. Such 'untrapping' cycles were performed for 
all databases until the low-energy region of the landscape con- 
verged. 

The methods described above are implemented in the pro- 
gram PATHS AMP LE, 38 which is a driver for OP T IM^as well 
as a tool for analysis of the resulting kinetic transition net- 
works. 
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2.4 Topological Characterisation 

To determine the topology of a given structure, it is first neces- 
sary to trace chains of head-to-tail dipoles within the cluster. 
An intuitive and robust method for identifying connectivity 
has been described by Miller and Wales. 14 Beginning with a 
particle /, the next particle in the chain is the one with the low- 
est (most favourable) dipole-dipole interaction energy with /, 
located in the half- space into which the dipole at / points. Sim- 
ilarly, the previous particle is the one with the lowest interac- 
tion energy in the complementary half-space. Connections in 
the chains that this procedure identifies will be referred to here 
as 'bonds.' 

In some structures, the bonds define a closed-loop chain that 
incorporates all the particles. If this loop cannot be unravelled 
into a trivial circle containing no crossings without breaking 
any bonds, then the structure is a knot. A knot is classified in 
the Rolf sen notation 39 according to the two-dimensional pro- 
jection of the chain that contains the smallest possible number 
of crossings. For the structures encountered in this work, the 
topology can be determined from an arbitrary projection by 
considering all combinations of splittings at the crossings and 
thereby evaluating the knot's Jones polynomial. 40 The Rolf- 
sen symbol can then be determined from tables.^ 
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Fig. 2 Upper panel: A path between the hexagonal antiprism and 
the distorted icosahedron for the Sti3 cluster at \i = 1.42. The graph 
shows the energy divided by the pair energy, V/e*, as a function of 
the integrated path length, s. Lower panel: Structures at labelled 
points on the path: top, in terms of the dipolar particles discussed in 
the text, and bottom, the dipole network. (The central particle, 
which does not participate in dipole bonding, is not shown.) 



3 Results and Discussion 



3.1 Sti3 clusters 

For Sti3, four structures successively become the global min- 
imum as fl is increased from zero: a distorted icosahedron, 
with D3 symmetry; a hexagonal antiprism, with Z)&/ symme- 
try; stacked rings of six and seven particles, with C s symme- 
try; and a planar thirteen-particle ring, with D\^h symmetry. 
The point groups here refer to symmetry operations that map 
the particle positions onto each other without considering the 
direction of the dipoles. 

Databases of minima and transition states have been con- 
structed for the three values of \l at which the global minimum 
changes, 1.42, 2.54, and 2.66, to show competition between 
the different morphologies, as well as at fl = 0, 0.2, 1, and 3.6 
to investigate the evolution of the potential energy landscape 
as a function of dipole strength. 

The corresponding disconnectivity graphs are presented in 
Fig.[T] For ji = 0, the expected single-funnel 'palm tree' mo- 
tijtETf s reproduced (Fig. [jj). This nomenclature refers to the 
fact that the branches of the graph are vertically rather short 
on the energy scale of the pairwise potential and that they are 
almost all directly connected to the central stem representing 
the basin of the global minimum. These features are indicative 
of multiple sequences of minima separated by relatively low 
barriers, converging on the lowest-energy structure. 



Introduction of a weak dipole causes a rapid increase in the 
number of minima on the landscape, owing to the multiplic- 
ity of stable arrangements of dipoles for each LJ morphol- 
ogy that are easily interconvertible via low barriers (Fig. [TJ).) 
Aside from the bifurcation of the branches that these arrange- 
ments cause, the form of the graph remains the same. When 
jl is increased to 1 , elements of this splitting are still present 
(Fig.[T]:). Two stable structures with practically identical ge- 
ometries and energies, but different arrangements of dipoles, 
lie substantially lower in energy than any other structure. Mor- 
phologies that do not exist on the LJ landscape have begun to 
emerge and, in terms of the optimal pair energy £*, the funnel 
appears to be shallower. 

At \i — 1.42, the global minimum changes from the dis- 
torted icosahedron to the hexagonal antiprism (Fig. [TJi.) The 
disconnectivity graph is still, broadly speaking, a palm tree 
but now has two low-lying minima that are separated by a 
substantial barrier. The lowest-energy single-step rearrange- 
ment in the database that connects the corresponding min- 
ima is presented in Fig. [2] The coordinate s in this plot is 
the integrated distance along the pathway between the initial 
and final structures in the 3A/^-dimensional Euclidean space of 
the particle positions. The two portions of a pathway from 
a transition state to its directly connected minima are defined 
by the approximately steepest-descent minimisations initiated 
from the transition state along the unstable mode. The ori- 
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Fig. 1 Upper panel: disconnectivity graphs for St 13 at ju = (a) 0; (b) 0.2; (c) 1; (d) 1.42; (e) 2.54; (f) 2.66; (g) 3.6. Branches to the 100 lowest 
energy minima are shown [except in (b) where 300 are shown for proper comparison with (a)]. Lower panel: global minima for Sti3. From 
left to right, the distorted icosahedron, the hexagonal antiprism, stacked rings, and the single ring. Particles are drawn as translucent spheres 
with arrows pointing along the dipole vector. Different chains are distinguished by their colour. Single chains are depicted with smoothly 
changing colour. 



entations of the dipole are not included in s and pure rota- 
tions of dipoles therefore do not contribute to the apparent 
length of the path. The interconversion mechanism is the well- 
known diamond-square-diamond rearrangement (DSD), pro- 
posed by Lipscomb, 43 which describes rearrangements of bo- 
ranes, 44 carboranes, 45 and metallaboranesJ^Ell It is also an 
important mechanism for LJ clusters.^ The pathway con- 
tains two consecutive DSD rearrangements, which succes- 
sively cleave both rings of the hexagonal antiprism at dia- 
metrically opposite points. To visualise this mechanism it is 
helpful to think of the hexagonal antiprism as a ring of twelve 
edge-sharing triangles. The DSD rearrangement rotates two 



opposite segments, each of four triangles, with respect to one 
another about a line connecting their centres. In the final 
structure (Fig. [2]l) the two segments now clasp each other in 
a relative orientation that is orthogonal to where they started 
(Fig. 0i). 

Dipolar 'bonding' adds a new dimension to the DSD re- 
arrangement, illustrated in Fig. [3^. Pairs of dipoles begin 
aligned along opposite edges of the diamond face. As the 
long diagonal contracts and the short diagonal elongates, the 
approaching dipoles reorient so that in the product they are 
aligned along the new short diagonal. The effect is to 'break' 
the bonds along the edges and create a bond along the diag- 
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Fig. 3 Rearrangement mechanisms for St 13 clusters as 
minimum- transition state-minimum triples: (a) the 
diamond-square-diamond (DSD) rearrangement; (b) the 
butterfly-tetrahedron (BTj) rearrangement; (c) the closed 
butterfly-double tetrahdron (B C T^) rearrangement; (d) the double 
BT^ (BTjT^/B) rearrangement. Transparent tubes show the 
underlying geometry, and filled tubes show dipole bonds. 



onal of the diamond. In the rearrangement between the dis- 
torted icosahedron and hexagonal antiprism of Sti3, this pro- 
cess breaks one bond in each ring and forms an interchain con- 
tact between two of the free ends. Along the pathway there 
exists a single twelve-particle helix (Fig. [2]:) wrapped around 
the central particle on the downhill path to the distorted icosa- 
hedron, which has three closed loops of dipoles. 

At jl = 2.54, the global minimum changes again, from the 
hexagonal antiprism to two stacked rings (Fig.[lJ). Compared 
to the first change of global minimum at jl = 1.42, the land- 
scape at this second change has more double-funnel charac- 
ter, since each of the two low-lying minima has a number of 
other minima associated with its main branch. In addition, the 
thirteen-membered planar ring, which is the global minimum 
in the strong dipole limit, now appears as a metastable struc- 
ture connected to the main stem of the graph at high energy. 
The associated region of the energy landscape may be consid- 
ered as a third funnel. 

The pathway connecting the lowest two minima is more 
complex than at the previous change, containing five transi- 
tion states, and traversing a variety of links and coils. Starting 
from the stacked rings (structure III in Fig. [I]), the first three 
steps convert the unlink to a coil, the coil to a link, and fi- 
nally the link to a different coil. Each of these transformations 



occurs by the butterfly-tetrahedron rearrangement (BTj), pro- 
posed by Wales et al.^ It is somewhat similar to Johnson's 
edge-bridging mechanism 50 in terms of the centres of mass: a 
butterfly moiety closes to form a tetrahedron. However, in the 
context of dipolar particles, this underlying structural change 
allows recombination of chains, interconverting topological 
isomers. This process is illustrated in Fig. [3]). Similarly to 
the DSD rearrangement, dipoles begin aligned along opposite 
edges of the butterfly. In contrast to the DSD process, the long 
diagonal contracts, but the short diagonal does not elongate 
significantly. As the incoming particles approach, the dipoles 
reorient to be along the interparticle vector, and a bond forms 
along the nascent tetrahedral edge. Simultaneously the dipoles 
at the hinge vertices align along the opposite edge. In this 
fashion, parallel chain segments are converted into perpendic- 
ular ones, creating a crossing in the chain of dipoles while 
keeping the number of dipole bonds conserved. 

The remaining steps of the interconversion from stacked 
rings to hexagonal antiprism involve rearrangements related 
to BTj. First, a coil is converted to a structure containing a 
five-membered ring and a seven-membered ring with a parti- 
cle in the centre (5,1,7). This process is difficult to describe 
in purely geometrical terms, but insight may be gained by 
considering which bonds are formed or broken in the pro- 
cess. This analysis allows us to describe what we will call the 
closed butterfly-double tetrahedron rearrangement (B C T^), il- 
lustrated in Fig. [3J;. We begin with a face-capped or 'closed' 
butterfly arrangement of particles, where bonds are formed 
along the same edges as are found in the BTj process. The 
butterfly closes around this particle, one wing at a time, form- 
ing bonds between the wing-tips and the capping particle, and 
along the hinge. The result can be thought of as two face- 
sharing tetrahedra. Comparison of Fig. |3J: and(3j) reveals the 
similarity in the overall changes to the dipole network. In this 
way a particle is released from two face- sharing tetrahedra in 
the coil, which becomes the central particle in the 5,1,7 ring 
system. 

The final step, in which the 5,1,7 system is converted to the 
hexagonal antiprism, is a concerted double BTj, or BTjTjB, 
rearrangement, illustrated in Fig. [3ji. In terms of the parti- 
cle centres, a butterfly and a tetrahedron sharing the hinge 
edge simultaneously close and open, interconverting. It differs 
from the other examples in that the butterfly's dipole bonds 
are not opposite each other, but adjacent, and the tetrahedra 
have bonds only along the edge that is broken in the butter- 
fly. In the context of the SI43 structure, the rearrangement of 
this fragment effects a particle exchange between the rings. At 
the transition state of the whole structure, none of the tetrahe- 
dral bonds and one butterfly bond for each of the rings exist, 
forming a twelve-particle helix. 

The global minimum in the high-dipole limit is a single 
thirteen-particle ring. This structure first replaces the stacked 
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Fig. 4 Upper panel: a path between the ring and stacked rings 
minima of the Sti3 cluster at jj, = 2.66. Lower panel: structures of 
minima and transition states along the path. 



rings as the Stockmayer global minimum at jl = 2.66. The dis- 
connectivity graph (Fig.[lJ.) shows distinct funnels for the two 
competing structures, but the stacked ring funnel is associated 
with a larger region of configuration space. As jl increases, 
and the ring funnel becomes increasingly favourable, the re- 
gion of configuration space associated with the stacked rings 
is expected to serve as a kinetic trap to structural relaxation. 

A pathway containing three transition states, which con- 
nects these competing minima, is presented in Fig. [4] The 
first two steps flatten the ring into an ellipse and twist it about 
its centre to form a coil. Although much of the bending energy 
penalty is removed, and favourable contacts are made between 
particles in the adjacent chains, the sharp bends at the ends of 
the ellipse are sufficient to increase its energy with respect to 
the ring. As the ellipse twists and folds into a coil, it is locked 
in place by a mechanism reminiscent of a DSD rearrangement, 
which is illustrated in Fig. [5^. The final step in the path is the 
now familiar BTj rearrangement, which interconverts the coil 
and the stacked rings, removing a crossing from the chain. It 
is interesting to note that in the previous pathway, at a smaller 
value of /I, all of the BTj steps lowered the energy upon form- 
ing the more compact tetrahedral site. In the current pathway, 
the situation is reversed, reflecting the increasing stability of 
planar moieties as the dipolar contribution to the energy in- 
creases. 

At large /I, the ring is firmly established as the global min- 
imum. It is separated from the next lowest minimum, the el- 
lipse, by 0.59£*, and from the third lowest energy structure by 
0.78e*. 

3.2 St2i clusters 

The most remarkable structures that appear as global min- 
ima for Stockmayer clusters are knotted geometries. 14 The 
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Fig. 5 Chain-based mechanisms: (a) a diamond-square-diamond 
rearrangement; (b) a budding rearrangement; (c) a DSD 
rearrangement that facilitates chains moving past one another; 



smallest cluster where a knot becomes energetically optimal 
is St2i. The global minimum of this cluster is topologically a 
trefoil knot over the range 1.7 < jl < 2.9, beyond which a pair 
of stacked twelve and thirteen particle rings becomes more 
favourable. The disconnectivity graph at jl = 2.9 is presented 
in Fig. [6| and exhibits a pronounced double funnel. As with 
the graph at jl = 2.66 for Sti3, the wider funnel is associated 
with the more compact structure — in this case, the knot. 

A low-energy pathway connecting these structures is pre- 
sented in Fig. [7] The pathway contains eleven transition states, 
and mostly involves rearrangement of, and particle exchange 
between, rings in a link. The first such conversion, from an 
unlink with ten and eleven particle rings to a link with nine 
and twelve particle rings, is achieved in three transition states. 
First, the unlink rearranges to a coil by the BTj mechanism. 
This coil contracts to a more compact bundle in a similar fash- 
ion to the Sti3 ring. The coil itself then twists, facilitated by 
coupled inter-chain DSD rearrangements, which enable chains 
to move past one another. Finally, the coil converts to the 
nine- and twelve-particle link by a second BTj rearrangement, 
which conserves a twist in the larger ring (Fig.|7J;). 

Two further exchanges, to form links of eight- and thirteen- 
particle rings and then the link of seven- and fourteen-particle 
rings, occur by a new mechanism, whereby a particle is 
smoothly transferred between adjacent chains. This 'budding' 
rearrangement is illustrated in Fig. [5]). A three-particle sec- 
tion of a chain bends sharply, expelling a particle and forming 
a new connection between the free ends. The released particle 
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Fig. 6 Upper panel: disconnectivity graph for St2i at ji =2.90. 
Branches to the 1000 lowest energy minima are shown. Lower 
panel: structures of the stacked ring (I) and trefoil knot (II) minima. 



-26.0 



-26.2 - 



-26.4 



-26.6 






10.0 


20.0 

s/<7 


30.0 


40.0 


(b) 




(c) 


(d) 


(e) 




Fig. 7 Upper panel: a pathway between the stacked rings and trefoil 
knot minima of St2i- The figure is shaded according to the number 
of chains in the structure (light=l, dark=2.) Lower panel: selected 
structures from the path: (a) the stacked rings; (b) a coil; (c) a link (9 
and 12 particles); (d) a link (7 and 14 particles); (e) the trefoil knot. 



is incorporated into an adjacent chain, which opens as the first 
one closes. At the transition state, the five sites are arranged 



in a closed butterfly geometry. The process can be thought 
of as an edge-bridging particle passing to the opposite edge 
via a face-bridged transition state. The latter exchanges are 
separated by coiling of the larger ring around the smaller by 
means of DSD rearrangements, shown in Fig. Bfc, and the rings 
are disposed such that a final BTj rearrangement produces the 
trefoil knot. 



3.3 St38 clusters 

The LJ38 cluster is unusual in that the global minimum is a 
truncated octahedron, and not based on an icosahedron.EUThe 
landscape is double-funneled, 52 and the lowest energy icosa- 
hedral minimum, while significantly higher in energy than the 
global minimum, is associated with a much larger region of 
configuration space. As such, the system provides a useful 
test of global optimisation procedures. 

St38 exhibits a more complex knotted global minimum than 
the St2i trefoil in the range 1 .6 < fl < 2.2lSlThis St 38 knot has 
eight crossings in its minimal projection and has the topology 
819 in Rolfsen's notation.^ 

The Lennard- Jones truncated octahedron ceases to be the 
global minimum at jl = 0.8. For values of jl between 0.8 and 
1.6, the truncated octahedron and knot — both metastable — are 
close in energy. Databases of minima and transition states 
have been constructed at jl = 1.22, where the two morpholo- 
gies are closest in energy. The disconnectivity graph is pre- 
sented in Fig. [8] It has the form of a gently sloping single fun- 
nel, with a form somewhere between a 'weeping willow' and 
a palm tree. 53 (The willow tree has longer vertical branches 
connected to the main stem, indicating multiple descending 
pathways towards the global minimum, like the palm tree, but 
with higher barriers between adjacent minima.) The global 
minimum has D^ symmetry and is depicted in Fig. [8] 

The truncated octahedron and knot structures can be inter- 
converted by a pathway containing just three transition states. 
Multiple concerted surface rearrangements occur in the first 
two steps, an example of which is presented in Fig. [9^. The 
square faces of the truncated octahedron minimum permit the 
previously characterised 54 square-diamond-diamond-square 
(SDDS) mechanism to occur. The final step is an intershell 
transfer of a particle from the surface to the core, expanding 
the octahedron to a pentagonal bipyramid (Fig. [9J3). Since 
we are considering a relatively weak dipole, the absence of 
dipole-dominated mechanisms, such as BTj, is not surprising. 
The 819 knot is still compact, but dipole bonding stabilises the 
open faces required to arrange the particles in such a way that 
the dipoles may orient to form a knot. 
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Fig. 8 Upper panel: disconnectivity graph for St38 at jd = 1 .22. 
Branches to the 1000 lowest energy minima are shown. Lower 
panel: the D^ global minimum (I) and the O^ fee minimum (II), 
showing the core-shell structure, and the 819 knot minimum, 
showing the dipole chain. 



4 Conclusions 

Clusters of spherical particles bound by simple isotropic po- 
tentials like Lennard- Jones and Morse are generally dom- 
inated by compact — often icosahedral — structures. The 
anisotropic Stockmayer potential differs from such cases be- 
cause of the particles' tendency to form chains. The fact that 
chain-like motifs can readily be identified introduces the new 
consideration of topology when characterising structure and 
rearrangements. Since the connectivity of particles within the 
chains is not fixed (as it is in a polymer), rearrangements may 
alter the topology. 

Some rearrangement mechanisms in Stockmayer clusters 
resemble those found in clusters bound by other potentials. 
For example, the diamond-square-diamond rearrangement, 
which is common in clusters bound by isotropic potentials, 





Fig. 9 Steps in the rearrangement of St38 from a truncated 
octahedron (structure II in Fig. [8} to an 819 knot (structure III), (a) A 
square-diamond-diamond- square (SDDS) rearrangement on the 
surface of the truncated octahedron (structure II in Fig. [8}. (b) A 
core- shell rearrangement in the final step of the pathway. Filled 
tubes highlight the geometry of the rearranging sites. 



also occurs in Stockmayer clusters. However, this famil- 
iar mechanism gains additional significance because of the 
change in orientation of the dipole vectors, which can allow 
chains to be disconnected in one direction and to become con- 
nected in another. Similarly, edge-bridging mechanisms seen 
in other clusters can change the topology of Stockmayer clus- 
ters. 

In addition to mechanisms that change the connectivity be- 
tween chains while preserving the association of each particle 
with its chain, there are budding processes. Here, a particle 
is smoothly ejected from one chain and absorbed into another, 
thereby allowing chains to change in size. 

These mechanisms vary in importance with the strength of 
the dipole moment. A crucial consideration is whether or 
not the pathway conserves the number of head-to-tail connec- 
tions between dipoles. The diamond-square-diamond mecha- 
nism exchanges two parallel connections for one in an orthog- 
onal orientation, and so is more favourable at lower dipole 
strengths, where the loss of a dipole-dipole bond does not in- 
cur a large energetic penalty. The butterfly-tetrahedron mech- 
anism and some of its variants preserve the number of bonds 
while interconverting topological isomers, and become more 
important at higher dipole strengths. 

Stockmayer clusters are frustrated systems in which the 
isotropic Lennard- Jones part of the potential drives the struc- 
ture towards compact highly-coordinated arrangements, while 
the dipolar interactions favour chain-like motifs. The analy- 
sis of structures and rearrangement mechanisms that we have 
presented here for Sti3, St2i, and St38 illustrates a variety of 
topology-changing events that offer a compromise between 
the competing terms in the potential. 
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Experimentally, clusters of colloids with a dipolar compo- 
nent to their interactions are being studied with increasing 
levels of control, and are being characterised in greater de- 
tail. While chain- and loop-based structures currently predom- 
inate, it should be possible to explore more compact clusters 
as well. The knotted topologies and other interesting struc- 
tures that arise for Stockmayer particles are obtained by ad- 
justing the relative strengths of the dipolar and isotropic parts 
of the potential. This flexibility can be achieved in suspensions 
of dipolar colloids by introducing tunable, depletion-induced, 
isotropic attraction via an additional component in the suspen- 



sion.^ The isotropic attraction can also be controlled through 
the surface chemistry of the colloids. n In both cases, the re- 
pulsive core of the particles is likely to be considerably harder 
than the r -12 term in Eq. |l]). However, it is the presence of 
an isotropic attractive tail and its competition with the dipolar 
interactions that are crucial in generating a rich selection of 
topologies. 

A number of open questions remain in the context of clus- 
ters and molecules with non-trivial topology. These issues in- 
clude the maximum complexity (number of crossings) that can 
be attained in a given system, and the precise factors determin- 
ing the favoured topology where a choice exists. The influ- 
ence of the interparticle potential and any fixed connectivity 
between particles will be important in this regard, and these 
considerations make contact with other knot- forming systems, 
such as proteins and synthetic organic molecules. The study 
of idealised model systems, such as the work presented here, 
should be a fruitful way of of gaining insight into the general 
phenomenon of knot formation. 
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